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Abstract 



A study of strange particle production in i/^ charged current interactions has 
been performed using the data from the NOMAD experiment. Yields of neutral 
strange particles (Kg, A, A) have been measured. Mean multiplicities are reported 
as a function of the event kinematic variables Ei,^ and as well as of the 
variables describing particle behaviour within a hadronic jet: xp, z and p\. Decays 
of resonances and heavy hyperons with identified Kg and A in the final state have 
been analyzed. Clear signals corresponding to K*^, S*^, E!~ and TP have been 
observed. 

Key words: neutrino interactions, strange particle production 



^ Now at Univ. of Perugia and INFN, Perugia, Italy 
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1 INTRODUCTION 



The production of strange particles in neutrino interactions can provide a 
testing ground for the quark-parton as well as for hadronization models. Neu- 
tral strange particles can be reliably identified using the y°-likc signature of 
their decays (Kg 7r"'"7r~, A — > pvr" and A pvT"*") in contrast to most other 
hadrons which require particle identification hardware. It is noteworthy that 
all previous investigations of strange particle production by neutrinos have 
come from bubble chamber experiments [1-25] . No other technique has so far 
yielded results on this subject. However, previous bubble chamber experiments 
with (anti)neutrino beams suffered from the low statistics of their samples. 

The NOMAD experiment [26] has collected a large number of neutrino in- 
teractions with a reconstruction quality similar to that of bubble chamber 
experiments. The order of magnitude increase in statistics can be used to im- 
prove our knowledge of strange particle production in neutrino interactions. In 
this paper we present measurements of the yields of neutral strange particles 
(Kg, A and A), as well as the yields of K*^, E*^, S~ and E° in i/^ charged 
current (CC) interactions. These results are compared to the predictions of 
the NOMAD Monte Carlo (MC) simulation program. 

The results of the present analysis are the measurements of: 

(1) the production properties of neutral strange particles in z/^ CC interac- 
tions. This study will allow tuning of the parameters of the MC simula- 
tion programs in order to correctly reproduce the production of strange 
particles by neutrinos; 

(2) the contribution of strange resonances and heavy hyperons to the total 
number of observed Kg , A and A. This will allow a quantitative theoretical 
interpretation of the A and A polarization measurements in z/^ CC deep 
inelastic scattering (DIS) reported in our previous articles [27,28]. 
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2 EXPERIMENTAL PROCEDURE 



2.1 The NOMAD experiment 



The main goal of the NOMAD experiment [26] was the search for 
oscillations in a wide-band neutrino beam from the CERN SPS. The main 
characteristics of the beam are given in Table 1. This search uses kinematic 
criteria to identify CC interactions [29] and requires a very good quality 
of event reconstruction similar to that of bubble chamber experiments. This 
has indeed been achieved by the NOMAD detector, and, moreover, the large 
data sample collected during four years of data taking (1995-1998) allows for 
a detailed study of neutrino interactions. The full data sample from the NO- 
MAD experiment corresponding to about 1.3 million z/^ CC interactions in 
the detector fiducial volume is used in the present analysis. The data are com- 
pared to the results of a Monte Carlo simulation based on modified versions 
of LEPTO 6.1 [30] and JETSET 7.4 [31] generators for neutrino interactions 
(with and W'^ cutoff parameters removed) and on a GEANT [32] based 
program for the detector response. Strange particle production is described 
by the set of default parameters in JETSET. To define the parton content 
of the nucleon for the cross-section calculation we have used the GRV-HO 
parametrization [33] of the parton density functions available in PDFLIB [34]. 
The above description of the MC will be referred to as the default MC. For the 
analysis reported below we used a MC sample consisting of about 3 million 
events. 

Table 1 

The CERN SPS neutrino beam composition at the position of the NOMAD detector 
(as predicted by the beam simulation program [35]). 



Neutrino 
flavours 


Flux 


CC interactions in NOMAD 


{E,) [GeV] 


rel.abund. 


{E,) [GeV] 


rel.abund. 




24.2 


1 


45.3 


1 




18.5 


0.0637 


40.9 


0.0244 




36.6 


0.0102 


57.1 


0.0153 




28.7 


0.0025 


49.9 


0.0015 
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Fig. 1. A sideview of the NOMAD detector. 
2.2 The NOMAD detector 



For a study of strange particle production the tracking capabihties of a detec- 
tor are of paramount importance. The NOMAD detector (Fig. 1) is especially 
well suited for this. It consists of an active target of 44 drift chambers, with a 
total fiducial mass of 2.7 tons, located in a 0.4 Tesla dipole magnetic field. The 
drift chambers (DC) [36] , made of low Z material (mainly carbon) serve the 
dual role of a nearly isoscalar target for neutrino interactions and of the track- 
ing medium. The average density of the drift chamber volume is 0.1 g/cm^, 
very close to that of liquid hydrogen. These drift chambers provide an over- 
all efficiency for charged particle reconstruction of better than 95% and a 
momentum resolution which can be parametrized as 

dp 0.05 0.008 -p 

where the track length L is in metres and the track momentum p is in GeV/c. 
This amounts to a resolution of approximately 3.5% in the momentum range 
of interest (less than 10 GeV/c). Reconstructed tracks are used to determine 
the event topology (the assignment of tracks to vertices), to reconstruct the 
vertex position and the track parameters at each vertex and, finally, to identify 
the vertex type (primary, secondary, V^, etc.). A transition radiation detec- 
tor [37] is used for electron identification. The pion rejection achieved for 
isolated tracks is 10^ with a 90% electron identification efficiency. A lead-glass 
electromagnetic calorimeter [38] located downstream of the tracking region 

provides an energy resolution of 3.2%/wi?[GeV] © 1% for electromagnetic 
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showers and is essential to measure the total energy flow in neutrino interac- 
tions. In addition, an iron absorber and a set of muon chambers located after 
the electromagnetic calorimeter are used for muon identification, providing a 
muon detection efficiency of 97% for momenta greater than 5 GeV/c. 

The large statistics of the data combined with the good quality of event recon- 
struction in the NOMAD detector allows a detailed study of strange particle 
production in neutrino interactions to be performed. 

2.3 Event selection and identification procedure 

The NOMAD experiment has collected 1.3 x 10^ i/^ CC events and has observed 
an unprecedented number of neutral strange particle decays. Such a decay ap- 
pears in the detector as a y°-like vertex: two tracks of opposite charge emerg- 
ing from a common vertex separated from the primary neutrino interaction 
vertex (Fig. 2). The V"°-like signature is expected also for photon conversions. 




Fig. 2. A reconstructed data event containing 3 vertices identified as Kg decays 
by our identification procedure. The scale on this plot is given by the size of the 
vertex boxes (3x3 cm^). 



The selection procedure for the CC event sample used in this analysis has 
been described in [27] . 

Since the NOMAD detector is unable to distinguish (anti)protons from pi- 
ous in the momentum range relevant to this analysis, our identification 
procedure relies on the kinematic properties of a decay. 

For the identification a kinematic fit method has been used as described 
in [27,39]. This fit has been performed for three decay hypotheses: Kg 7r"'"7r~, 
A —>■ pvr", A —>■ pn'^ and for the hypothesis of a photon conversion to e+e". 
The output of the kinematic fits apphed to a given y° vertex consists of four 
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values of Xyo describing the goodness of these fits. Different regions in the 
four-dimensional Xyo space populated by particles identified as Kg, A and A 
have been selected. Identified V'^ are of two types: 

• uniquely identified V'^, which, in the four-dimensional Xyo space described 
above, populate regions in which decays of only a single particle type are 

present; 

• ambiguously identified which populate regions in which decays of differ- 
ent particle types are present. 

The treatment of ambiguities aims at selecting a given decay with the high- 
est efficiency and the lowest background contamination from other V° types. 
An optimum compromise between high statistics of the identified y° sample 
and well understood background contamination is the aim of our identifica- 
tion strategy. The MC simulation program has been used to define the criteria 
for the kinematic selection and to determine the purity of the final Kg, A 
and A samples. We selected a sample consisting of more than 90% of uniquely 
identified V^. Results are reported in Table 2. 

Table 2 

Efficiency (e) and purity (P) for each selected category. Numbers of identified 
neutral strange particles in the data are also shown in the last column. 





6(%) 


P(%) 


Data 




22.1 ±0.1 


97.2 ±0.1 


15074 


A 


16.4 ±0.1 


95.9 ±0.1 


8087 


A 


18.6 ±0.5 


89.7 ±0.7 


649 



The total sample in our data contains 15074 identified Kg, 8087 identified A 
and 649 identified A decays, representing significantly larger numbers than in 
previous (anti)neutrino experiments performed with bubble chambers [1-25]. 

Fig. 3 shows the invariant mass and cr distributions for identified Kg, A and 
A. The measured mass and the lifetime of identified neutral strange particles 
are in agreement with the world averages [40]. The corresponding results are 
given in Tables 3 and 4. 

The efficiencies and purities reported in Table 2 are momentum dependent. 
However, we have checked that they are applicable to the data because the mo- 
menta distributions of identified and of their decay products are identical 
in the data and MC simulation. 

In the rest of this paper we will always present efficiency corrected distribu- 
tions. 
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Table 3 



Measured mass and resolution (in MeV/cP) for both MC and data. 





MC 


Data 


PDG [40] 


Mass 


Resolution 


Mass 


Resolution 


Mass 




497.9 ± 0.1 


11.0 ± 0.1 


497.9 ± 0.1 


11.3 ± 0.1 


497.672 ± 0.031 


A 


1115.7 lb 0.02 


4.0 ± 0.03 


1116.0 it 0.05 


4.45 lb 0.05 


1115.683 it 0.006 


A 


1116.0 ± 0.1 


4.1 ± 0.1 


1116.3 ± 0.1 


4.8 ± 0.2 


1115.683 ± 0.006 



Table 4 

Measured ct (in cm) for a given type for both MC and data. 





MC 


Data 


PDG [40] 


s 


2.60 ± 0.01 


2.72 ± 0.03 


2.6786 


A 


7.91 ± 0.02 


8.07 ± 0.12 


7.89 


A 


7.82 ± 0.06 


7.33 ± 0.33 


7.89 
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Fig. 3. Invariant mass distributions and efficiency corrected ct distributions for 
identified Kg (top), A (center) and A (bottom) in CC DIS events for both data 
(points with error bars) and MC (histogram). 
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3 YIELDS OF NEUTRAL STRANGE PARTICLES 



We have studied the production rates of neutral strange particles (Kg, A, A) 
in Pfj^ CC interactions. The particles can be produced at the primary vertex 
and also from secondary interactions of primary particles with the detector 
material. Neutral strange particles produced via resonance or heavier hyperon 
decays are classified as primary V'^. We applied a correction obtained from 
the Monte Carlo to the yields of neutral strange particles in the data in order 
to extract the yields at the primary vertex. 



3. 1 Integral yields of neutral strange particles 



The measured yield per CC interaction for each type is defined as: 

where Nyo is the number of reconstructed and identified in the number 
Nvy,cc of reconstructed CC events and ^ is a correction factor calculated 
as: 

^ eyo X Br{V^ ^ h+h-)' 

where eu^cc = (85.30 ± 0.02)% is the efficiency (reconstruction and identi- 
fication) of CC events in the fiducial volume, and eyo is the global 
efficiency which takes into account the contribution from particles produced 
in the fiducial volume, but decaying outside. Pyo is the purity of the final 
sample, and Br{V^ h^h~) is the branching ratio for a given type 
decaying to a pair of charged hadrons. 

Table 5 

Integral yields of primary in CC interactions in both the data and in the 
default MC. The errors are only statistical. 



y° type 


^DATA 


(%) 


T^MC /rpDATA 




6.76 ± 0.06 


9.50 ± 0.02 


1.40 ±0.01 


A 


5.04 ± 0.06 


8.10 ±0.02 


1.61 ±0.02 


A 


0.37 ±0.02 


0.60 ±0.01 


1.62 ±0.03 



Table 5 shows the overall inclusive production rates for Kg, A and A in z/^ 
CC interactions compared to the MC predictions. Note that the production 
rates in the default MC are a factor of 1.4 to 1.6 higher than in the data. This 
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could be explained in part by the choice of the LEPTO [30] and JETSET [31] 
parameters in the NOMAD event generator. We have to emphasize that the 
so-called ss suppression factor (PARJ(2) parameter in JETSET) - defined as 
the ratio of the probability 7^ of producing an ss pair to the probability 7^ (7^) 
of producing a uu (dd) pair in the fragmentation chain - was set to its default 
value of 0.3 in the default NOMAD MC production[3]. However, this parameter 
was measured to be about 0.2 in previous bubble chamber experiments: for 
example, the values obtained by the BEBC WA21 Collaboration [41] in a 
neutrino beam similar to ours are 0.200 ± 0.022{stat) ± O.OlO(sys) for up and 
0.207 ± 0.018{stat) ± 0.020{sys) for up interactions. Moreover, later results 
from OPAL (0.245) [42], DELPHI (0.23) [43], E665 (0.2) [44], ZEUS [45] and 
HI [46] collaborations support a value close to 0.2 for this parameter. 

Table 6 

Numbers of u^ CC events with a specified combination of observed neutral strange 
particle decays for both the default MC and data. X indicates the hadronic system 
accompanying the observed . 



Channel 


Number of observed events 


mc/data 


MC 


Data 


AX 


11686 


7778 


1.50 ±0.02 




18971 


14228 


1.33 ±0.01 


AX 


831 


594 


1.40 ± 0.07 


KOROX 


485 


284 


1.7 ±0.1 


AK^X 


617 


247 


2.5 ±0.2 


AAX 


98 


40 


2.5 ±0.4 


K^AX 


24 


15 


1.6 ±0.5 


AAX 


19 


10 


1.9 ±0.7 


ak°kOx 


7 


2 


3.4 ±2.6 


K°KOKOX 


2 


4 


0.6 ±0.4 



However, the problem of the inaccurate description of neutral strange particle 
production in the MC is a more complex one. This is illustrated in Table 6, 
where we give the observed numbers of z/^ CC events for 10 exclusive multi-y° 
channels in the data compared to the default MC predictions. The number of 
MC events in Table 6 is renormalized to the same number of u^ CC events 
as in the data. From this comparison one can conclude that it is not possible 

^ The value 0.3 has been suggested by the authors of JETSET as the default for 
this parameter. 
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to rescale just a single parameter (the ss suppression factor) in order to de- 
scribe the neutral strange particle production observed in the data. Rather, 
the discrepancy is due to a combination of several parameters which describe 
the probability that an s(s)-quark appears as a meson/baryon(antibaryon), 
the probability that a strange meson/baryon(antibaryon) appears electrically 
neutral, etc. A tuning of the JETSET parameters to reproduce the yields of 
neutral strange particles observed in the NOMAD data is a subject of an 
analysis currently in progress. 

The integral yields reported in Table 5 can be compared to previous measure- 
ments summarized in Table 7. The Kg rates from Table 5 have been converted 
into K^{— + K^) rates by multiplying by a factor of 2. 

For completeness we have recalculated the integral overall yields taking into 
account contributions from both primary and secondary These results are 
given in Table 7 and denoted by a star (★). Our overall yields are consistent 
with the results of the i/-Ne experiment [19] performed in a similar neutrino 
beam. However, our primary yields of Kg and A are about 30% lower and the 
primary yield of A is ~20% lower. 



3.2 Yields of neutral strange particles as a function of kinematic variables 

To investigate neutral strange particle production mechanisms we have mea- 
sured the average Kg, A and A yields as a function of the neutrino energy E^, 
the invariant effective mass squared W"^ of the hadronic system, the invariant 
square of the four- momentum transfer from the neutrino to the target Q^, and 
the Bjorken scaling variable xbj- 

The yields of Kg , A and A are shown in Fig. 4. The A yield shows a behaviour 
which is almost independent of E^, W"^ and Q"^ after a sharp initial rise. It 
drops at large values of XBj- On the other hand, the yield of Kg rises steadily 
with Ey and W^^, reaches a plateau at large and falls with increasing XBj- 
Similar observations have been made by previous experiments, but with larger 
statistical uncertainties. The A yields as a function of kinematic variables are 
measured for the first time in a neutrino experiment. In general they show a 
behaviour similar to that of the Kg. However, as expected, clear W'^ and E^ 
thresholds are present in the A production. 

The ratios of yields for Kg /A and A/ A are presented in Fig. 5. 
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Table 7 

Inclusive yields of neutral strange particles in CC interactions measured in this 
analysis and in previous bubble chamber experiments. Nk, N\ and iV^ are the ob- 
served numbers ofK^, A and A, respectively. stands for + K^. See text for 
explanation of a star (*). 



Reaction 






rate 


Na 


A rate 




A rate 


fErfl 

L J 


(GcV) 




(Vi) 




(Vi) 




(Vi) 


NOMAD 


45 


15075 


1 3 52 + 1 2 


8087 




649 


f) 07 _|_ f) n9 
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ooy 
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± / .0 It u.y 


ioU 


/i c: 4- n /I 
4.0 It U.4 
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v - n [14] 


62 


234 


20.8 lb 1.6 


157 


7.1 lb 0.7 






u - p [14] 


62 


154 


17.7 lb 1.6 


77 


4.3 lb 0.6 






u - n [17] 


62 




20.5 lb 1.1 




6.6 lb 0.7 






u - p [17] 


62 




17.4 lb 1.2 




4.4 lb 0.5 






u - A [21] 


~10 


82 


7.1 ±0.8 


76 


3.1 lb 0.4 






u - p [10] 


~45 


23 


15 lb 4 










u - p [15] 


~50 








7.0 ± 0.8 






u - n [15] 


~50 








7.0 lb 1.2 







3.3 Comparison with LUND model predictions 



The measured yields of K^, A and A particles are compared in Figs. 6, 7, 8 
to the predictions of the default NOMAD MC simulation (see section 2.1). 
A reasonable agreement in shapes is observed, while the discrepancy in the 
overall normalization is about a factor of 1.3 to 1.5. 

The kinematic variables E^, W'^, and xsj are not independent. So, for 
example, discrepancies between the data and MC at high are reflected in 
discrepancies at low xbj- 
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Fig. 4. Measured yields of the Kg, A and A as a function of E^, W"^, and xbj- 
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Fig. 5. Ratios of measured yields for Kg /A and A/A as a function of E^, , 
and X Bj ■ 
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Fig. 6. Yields in the default MC (histogram) and in the data (points with error bars) 
forK^ as a function of Ei^, W"^, and xbj- The MC/Data ratios and their fit to 
a constant are also shown. 
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Fig. 7. Yields in the default MC (histogram) and in the data (points with error bars) 
for h. as a function of Ey, W^, and xbj- The MC/Data ratios and their fit to a 
constant are also shown. 
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Fig. 8. Yields in the default MC (histogram) and in the data (points with error bars) 
for A as a function of Ey, W"^ , and xbj- The MC/Data ratios and their fit to a 
constant are also shown. 
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4 PRODUCTION PROPERTIES 



We have also performed a detailed analysis of kinematic quantities describing 
the behaviour of neutral strange particles (Kg, A and A) inside the hadronic 
jet. This study allows an investigation of the dynamics of fragmentation. The 
differences in the production properties of Kg, A and A are seen most clearly 
here. The following distributions are of interest: xp = 2p*^/W (Feynman-x is 
the longitudinal momentum fraction in the hadronic center of mass system), 
the transverse momentum squared, of a particle with respect to the current 
(hadronic jet) direction and the fraction z = Eiab{V^)/Eiah{si\\ hadrons) of 
the total hadronic energy carried away by the neutral strange particle in the 
laboratory system. 



4-1 xp distributions 



■g 0.25 - 
I 




Fig. 9. Efficiency corrected xp distribu- Fig. 10. Efficiency corrected xp distribu- 
tion for Kg mesons. tion for A hyperons. 

The efficiency corrected xp distributions observed in the data for neutral 
strange particles are shown in Figs. 9, 10, 11. These distributions indicate 
that A are produced mainly in the target fragmentation region {xp < 0), 
while Kg are peaked in the central region with an asymmetry in the forward 
direction. A are produced in the central xp region {\xp\ < 0.5). One way 
to quantify differences in the xp distributions is to define an asymmetry pa- 
rameter A — (Np — Nb)/{Np + Nb), where Np and Nb are the numbers of 
particles produced forward and backward, respectively, in the hadronic centre 
of mass. The asymmetry parameters A and mean values of xp in both data 
and MC are listed in Table 8. They are consistent with previous measure- 
ments [11,12,14,19,23-25]. The observed disagreement between data and MC 
is probably due to the fact that the MC simulation does not properly describe 
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Fig. 11. Efficiency corrected xp distribution for A hyperons. 
the relative contributions of different V'^ production mechanisms. 
Table 8 

Mean values (xp) and asymmetry parameters A of the xp distributions for Kg, A, 
A in both MC and data. 





MC 


DATA 


{xf) 


A 


{xf) 


A 




0.055 ± 0.001 


0.152 ± 0.002 


0.064 ± 0.001 


0.256 ± 0.004 


A 


-0.296 ± 0.001 


-0.649 ± 0.002 


-0.295 ± 0.002 


-0.589 ± 0.004 


A 


0.006 ± 0.002 


-0.03 ± 0.01 


0.04 ± 0.004 


0.18 ±0.02 



^.2 z distributions 



Efficiency corrected z distributions for Kg, A and A are shown in Figs. 12, 
13, 14. A turn-over at small values of z can be seen for A hyperons, but not 
for Kg and A. A turn-over at small values of z for Kg and A was observed in 
some of the previous neutrino experiments [8,11,14,21] and was not observed 
in others [10,12]. We note that uncorrected z distributions show such a turn- 
over for all y° types in our experiment as well, due to a less efficient 
reconstruction at low momenta. 

Below we study separately z distributions of neutral strange particles produced 
aX Xf < Q (see § 4.2.1) and at xj? > (see § 4.2.2). Detailed information on 
the mean values of z distributions for both data and MC is given in Table 9. 
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Fig. 12. Efficiency corrected z distribu- Fig. 13. Efficiency corrected z distribu- 



tion for Kg . 
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Fig. 14. Efficiency corrected z distributions for A. 



4-2.1 z distributions in the target fragmentation region 



The efficiency corrected z distributions of Kg, A and A measured in tfie tar- 
get fragmentation region are sliown in Fig. 15 (left). One can see tliat tlie z 
distribution of Kg mesons has a maximum at 2; — > and decreases faster at 
larger z (compared with Fig. 12) and that the Kg mesons produced in the 
target fragmentation region carry in general a small fraction of the hadronic 
jet energy. A hyperons are believed to be produced mostly from the remnant 
di-quark fragmentation and the shape of the z distribution is similar to that 
shown in Fig. 13. The turn-over in the z distribution is observed for A hyperons 
produced in the target fragmentation region. 
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Table 9 

Mean values of z distributions for Kg, A, A measured for the full sample and for 
xf <0 and xp > regions in both MC and data. 







full sample 


XF <0 


XF > 




MC 
DATA 


0.218 ±0.001 
0.226 ± 0.001 


0.092 ± 0.001 
0.105 ±0.001 


0.312 ± 0.001 
0.299 ±0.001 


A 


MC 


0.227 ± 0.001 


0.179 ± 0.001 


0.462 ± 0.003 




da:i_a. 


0.250 ± 0.002 


0.206 ± 0.001 


0. 13-1 ± 0.005 


A 


MC 


0.215 ± 0.003 


0.139 ± 0.002 


0.296 ± 0.004 


data 


0.242 ± 0.005 


0.147 ±0.005 


0.308 ± 0.008 



4-2.2 z distributions in the current fragmentation region 

The efficiency corrected z distributions of , A and A measured in the current 
fragmentation region are shown in Fig. 15 (right). This kinematical region 
is interesting because of the it or d (anti) quark fragmentation into A or A 
hyperons. All three z distributions show similar behaviour but with different 
mean values of z. The z distribution of A hyperons at a: i? > is drastically 
different from that in the target fragmentation region. This is evidence for 
the fragmentation of the outgoing u quark into a A hyperon. In fact, the z 
distribution of A in the current fragmentation region is a measure of the D^{z) 
fragmentation function (normalized to unity in Fig. 15). The z distribution of 
A hyperons is sensitive to the D^{z) fragmentation function with a possible 
contribution from the D^{z) process. One can see that it is harder than the 
one measured in the target fragmentation region. 



4-3 Discussion 

There are different mechanisms responsible for Kg, A and A production in 
the neutrino CC DIS process which are expected to give different xf and z 
distributions for these particles. 

- The xp distribution of Kg mesons produced promptly in the W~^d — > u pro- 
cess that requires at least two quark- antiquark pairs to be created {dd and 
ss) is expected to be central. A contribution from heavier strange particle 
decays (mainly from K*~^) produced from the fragmentation of the outgo- 
ing u quark can result in a forward xp distribution for Kg mesons. Also Kg 
mesons from a fragmentation of the outgoing (anti)quark in W'^d c ^ s 
and W'^u — > d processes are expected to be produced in the forward xp 
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region and to carry a larger fraction of the jet energy. 

- A hyperons can be produced from the fragmentation of the micleon di-quark 
remnant promptly and via the decay of heavier strange baryons at xp < 0- 
A hyperons can be produced also ai xp > from the outgoing u quark 
fragmentation. 

- The production of A hyperons in neutrino scattering from a valence quark 
requires three quark- ant iquark pairs to be created {uu, dd and ss) and is 
expected to populate the central region of the xp distribution. There could 
also be a contribution from the outgoing antiquark fragmentation into a A 
hyperon (in the W^u d ^ A process) which can produce these baryons 
in the forward xp region. 

4.4 Pt distributions 

The efficiency corrected distributions of Kg, A and A in the data are pre- 
sented in Figs. 16, 17, 18. They show an exponential behaviour of the form 
Cexp{—B ■ p^) in the region p|> < 0.5 (GeV^/c^) and a deviation from this 
dependence at higher p^. We measured the slope parameter B in the region 
< Pt < 0.5 GeV^/c^ for each category in three kinematic regions: all xp, 
Xp <{i and Xp > 0. The results are hsted in Table 10. The values of the slope 
parameter found in different xp regions are similar, except for A hyperons. 

Table 10 

The slope parameter B{GeV/c)~'^ of the distribution for Kg, A, A measured 
separately for the full sample and for xp < Q and > regions in both MC and 
in data. 





MC 


DAm 




full sample 


< 


> 


full sample 


xp <0 


xp>0 


s 


5.72 ±0.03 


5.61 ±0.04 


5.79 ±0.03 


5.21 ±0.10 


5.40 ±0.21 


5.15 ±0.11 


A 


4.35 ±0.03 


4.30 ±0.03 


4.58 ±0.06 


4.12 ±0.13 


4.18 ±0.15 


4.07 ±0.27 


A 


3.89 ±0.10 


4.10 ±0.13 


3.70 ±0.14 


4.42 ± 0.47 


6.59 ± 0.74 


3.30 ± 0.64 



4-5 (pf.) versus xp distributions 

We have also studied the dependence of the average {p\) on xp for K^, A and 
A (see Figs. 19, 20, 21). For the first time in a neutrino experiment the good 
quality of event reconstruction combined with the large statistics of the data 
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collected allows the study of these distributions for neutral strange particles^. 
The observed discrepancy between the data and simulated events in the region 
xp > 0.3 could be attributed to the absence of QCD effects in our Monte Carlo 
simulation program: the so-called soft-gluon effect could change the leading 
particle [xp —>■ 1) behaviour inside the hadronic jet since the forward scattered 
quark is strongly accelerated and is therefore expected to radiate gluons, thus 
broadening the forward p'^ distribution. It has also been verified that there 
is no accumulation of misidentified in the region where the disagreement 
between the MC simulation and the data is observed. 



^ Similar distributions obtained for charged particles in bubble chamber neutrino 
experiments [47] have been used to tune Monte Carlo simulation programs. 
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Fig. 15. Efficiency corrected z distributions for xp < (left) and for xp > (right) 
/or Kg (top), A (centre) and A (bottom). 
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Fig. 16. Efficiency corrected distribu- Fig. 17. Efficiency corrected distribu- 
tion for K^. tion for A. 
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Fig. 18. Efficiency corrected p\ distribution for A. 
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Fig. 19. Efficiency corrected (p^) versus Fig. 20. Efficiency corrected {p^) ver- 
XF distribution for Kg in data (points susxp distribution for K in data (points 
with error bars) and MC (histogram). with error bars) and MC (histogram). 
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Fig. 21. Efficiency corrected {p^) versus xp distribution for A in data (points with 
error bars) and MC (histogram). 
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5 STRANGE RESONANCES AND HEAVIER HYPERONS 

Apart from many good physics reasons a study of the production of resonances 
and heavy hyperons is also of great importance for tuning the LUND model 
parameters and for the theoretical interpretation of the A and A polariza- 
tion measurements reported in our previous articles [27,28]. This is essential 
because A hyperons originating from the decays E* — > Att, — > A7 and 
S Att inherit a polarization from the parent particles and this polariza- 
tion is different from that of a directly produced A. Information about E°, S, 
E*. E* and K* yields can be obtained from an analysis of their decays into 
channels containing identified neutral strange particles [39,48]. 

Previous bubble chamber experiments with (anti) neutrino beams suffered from 

a lack of statistics. For example, the BEBC WA21 Collaboration [24] reported 
the observation of 149 ± 29 42 ± 19 K*' , 134 ± 19 E*+ and less than 

10 E*~ in i/fj^p CC interactions, while the Fermilab 15-ft bubble chamber E380 
Collaboration [19] found 94 ± 25 E° and 4 in i/^Ne CC events. 

5.1 A procedure for signal extraction 

Our aim was to extract the fraction of neutral strange particles which are 
decay products of resonances and heavier hyperons from the corresponding 
invariant mass distributions. To construct such distributions we combine the 
neutral strange particle with all possible charged tracks (of appropriate sign) 
emerging from the primary vertex except those identified as muons or elec- 
trons. We have also studied the (A 7) combinations, where photons are iden- 
tified as conversions in the detector fiducial volume via our identification 
procedure. The resulting distributions are fitted by a function describing both 
the combinatorial background and the resonance signal. 

The combinatorial background (BG) can be approximated by any function of 
the form: 



where Pnim — Mtu) is a polynomial of order n vanishing at m = Mth, Mtu ~ 
Myo -|- m^r, and Tail{m) is any function vanishing at m — > 00 faster than P„ 
increases. 

We have chosen the following BG paramctrization: 



Pn{m-Mth)-Tail{m), 



(2) 



(3) 



where A — m — Mi 



th- 
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For a resonance signal the standard relativistic Breit-Wigner (BW) func- 
tion [49] is used: 

BW{m)^ ^ , J,^ ,,,^J -], with T = To(^]^ (4) 

where Mq, Fq are the resonance mass and width, respectively, and q is the 
momentum of the decay product in the resonance rest frame {qo corresponds 
to Mo). 

Finally, we have fitted the invariant mass distributions by: 

^ = BG(A) + a,BW'(m), (5) 

for all combinations except (A tt^), where two peaks due to S*~ — > Att" and 
S~ —>■ Att" decays are expected. Here BW'{m) is the Breit-Wigner function 
of equation (4) normalized to unity. Such a fit is valid in all cases when the 
experimental mass resolution is small compared with the natural width of the 
resonance. 

Similarly, for the (A tt") case we have used 
dN 

— = BG(A) + a5BW^.- (m) + ueBWi- (m) (6) 
dm 

where the invariant mass resolution is used for the width Fq in the Breit- 
Wigner function corresponding to the S~ decay. In the above formulae ai to 
are parameters of the fit. 



In such an approach using the HESSE and MINOS procedures of MINUIT 
the parameter 05(06) gives the number of signal events with the corresponding 
error which takes into account possible correlations between different param- 
eters. 

As a consistency check we have also tried an alternative approach which was 
to fit the invariant mass distributions with: 

^^{l + a,BW{m))BG{A), (7) 
am 

and, similarly, with 
dN 

— = (1 + a^BWj,*- (m) + UgBWe- (m))5G'(A), (8) 
dm 

and extracted the corresponding number of signal events. 

The results obtained using these two approaches were found to be similar. In 
what follows we present our results using the first method. 
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5.2 Results 



The yields of resonances and heavy hyperons have been studied in different 
kinematic regions and for neutrino interactions on different target nucleons. 

In NOMAD it is to some extent possible to separate neutrino interactions on 
the neutrons and protons by imposing a cut on the sum of charges (Qtot) oi 
all the outgoing tracks at the primary neutrino interaction vertex. 

We select up events requiring Qtot > 1- According to the MC simulation, in this 
proton-like sample 76% of the events are true up interactions. The un events 
are selected by the requirement Qtot ^ 0. The purity of the corresponding 
neutron-like sample is about 85%. 

Since the hadron production mechanisms in the target and in the current 
fragmentation regions are expected to be different, it is important to study 
separately the yields of resonances and heavy hyperons at Xi;' < and xp > 0. 
Such a study is also necessary for a correct theoretical interpretation of the A 
(A) polarization measurements reported in our previous papers [27,28]. 

In the following we denote as MC(pred.) the true number of heavy strange par- 
ticles reconstructed in the MC, and MC(meas.) the number of heavy strange 
particles extracted from the MC sample using our fitting procedure. Both 
quantities are normalized to the number of z/^ CC events in the data. 

Note that MC(pred.) and MC(meas.) can be slightly different due to hmita- 
tions of the signal extraction procedure described in section 5.1. The threshold 

and smearing effects in the invariant mass distributions are at the origin of this 
discrepancy. The ratio MC(pred.)/MC(meas.) will therefore be used to correct 
the yields of heavy strange particles extracted from the data (see section 5.3). 

The following resonances and heavier hyperons have been studied in the 
present analysis. 



5.2.1 K*± 

The fitted (Kg tt^) invariant mass distributions are shown in Fig. 22 for both 
MC and data samples. Detailed information on the number of extracted K*^ 
events and the K*^/Kg ratio is given in Tables 11 and 12. For the K*^ mass 
and width we have used 891.66 MeV and 50.8 MeV respectively. The Qq value 
is 291 MeV/c [40]. 

It is interesting to note a more abundant K*^ than K*~ production in u^ 
CC DIS. This can be explained by the fact that the outgoing u quark can 
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Fig. 22. KgYr"*" (top) and Kg7r~ (bottom) invariant mass distributions for both MC 
(left) and data (right). The solid lines are the results of the fit, while the dotted lines 
describe the background term. In the MC plots the additional histograms refer to the 
reconstructed true heavy strange particles. 

fragment directly into a K*~^, while both u and s quarks needed to produce a 
K*~ meson have to be created in the fragmentation process. 

One can see that there is a significant difference between K*''' yields in the 
default MC simulation and the NOMAD data (about a factor of 2). 



5.2.2 E*± 

For the S*"*^ mass and width we have taken the values from [40]: m(S*''')=1382.8 
MeV, r(S*+)=35.8 MeV, m(S*-)=1387.2 MeV, r(E*-)=39.4 MeV. The qo 
value is 208 MeV/c. 
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Table 11 



K*^ Kg7r+ summary 



N(K*+) 


full sample 


Kg fragmentation region 


type of target nucleon 








> 


vp 


vn 


DATA 


2036 ± 121 


315 ± 56 


1731 ± 108 


1006 ± 87 


1032 ± 84 


MC(meas.) 


5373 ± 104 


726 ± 47 


4744 ± 93 


1963 ± 67 


3516 ± 80 


MC(pred.) 


5953 


886 


5067 


2206 


3748 


N(K'^+)/N(K0) 


uncorrected 


DATA (%) 


13.5 ±0.8 


9.7 ± 1.7 


14.6 ± 0.9 


15.7 ± 1.4 


11.9 ± 1.0 


MC (%) 


27.3 ±0.5 


14.5 ±0.9 


31.6 ±0.6 


29.0 ± 1.0 


26.5 ± 0.6 


Table 12 

K*^ — KJ^tt^ sinrnriarij 


N(K*-) 


full sample 


Kg fragmentation region 


type of target nucleon 






XF < 


XF > 


vp 


un 


DATA 


1146 ± 89 


288 ± 44 


865 ± 78 


377 ± 52 


775 ± 73 


MC(meas.) 


2304 ± 74 


639 ± 38 


1664 ± 63 


729 ± 39 


1576 ± 63 


MC(pred.) 


2467 


723 


1743 


734 


1733 


N(K*-)/N(KO) 


uncorrected 


DATA (%) 


7.6 ±0.6 


8.9 ± 1.3 


7.3 ±0.7 


5.9 ±0.8 


8.9 ±0.8 


MC (%) 


11.5 ±0.4 


12.7 ±0.8 


11.1 ±0.4 


10.8 ±0.6 


11.9 ±0.5 



The fitted invariant mass distributions for (A vr^) combinations in both MC 
and data samples are shown in Fig. 23. Detailed information on the number 
of extracted S*^ events and the E*^/A ratio is given in Tables 13 and 14. 



A striking difference between the S* yields in the default MC simulation and 
in the NOMAD data (a factor of about 3 or even larger) is observed. 

5.2.3 E- 

For the S~ mass we have used 1321.32 MeV [40], for the width the experi- 
mental resolution of 10 MeV has been taken. The go value is 139 MeV/c. 

The bottom plots in Fig. 23 show evidence for Avr" decays. Detailed 

information on the number of extracted S~ events and the S~/A ratio is given 
in Table 15. 
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Fig. 23. Avr"'' (top) and Att" (bottom) invariant mass distributions for both MC 
(left) and data (right). The solid lines are the results of the fit, while the dotted 
lines describe the background term. In the MC plots the additional histograms refer 
to the reconstructed true heavy strange particles. 



5.2.4 S° 



For the E° mass we have taken the value from [40]: m(E°)=1192.6 MeV, while 
for the width the experimental resolution of 9 MeV has been used. The go 
value is 74 MeV/c. The E° peak has been fitted by a Gaussian function. 

Fig. 24 shows the fitted invariant mass distributions for (A 7) combinations 
in both Monte Carlo and data samples. The corresponding photons have been 
reconstructed as conversions in the DC fiducial volume and identified by our 
identification procedure. The quality of the photon reconstruction is illus- 
trated by the (77) invariant mass distributions shown in Fig. 25 for both MC 
and data: a peak corresponding to the 7r° signal is evident. 
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Table 13 



S*"*" Avr^ summary 



N(E*+) 


full sample 


A fragmentation region 


type of target nucleon 






XF <0 


XF >0 


vp 


un 


r> ATA 


416 ± 80 


358 ± 65 


63 ±47 


297 ±61 


120 ±51 


MC(meas.) 


2070 ± 68 


1427 ± 57 


649 ± 37 


1321 ± 49 


754 ± 46 


MC(pred.) 


1783 


1254 


529 


1150 


634 


N(S'^+)/N(A) 


uncorrected 


DATA (%) 


5.2 ± 1.0 


6.4 ± 1.2 


2.5 ± 1.9 


8.6 ± 1.8 


2.6 ± 1.1 


MC (%) 


17.0 ±0.6 


15.9 ±0.6 


20.6 ± 1.2 


32.6 ± 1.2 


9.3 ±0.6 



Table 14 



S* Avr summary 



N(E*-) 


full sample 


A fragmentation region 


type of target nucleon 






XF <0 


XF > 


vp 


vn 


DATA 


206 ± 63 


121 ± 51 


93 ±37 


100 ± 35 


111 ±52 


MC(meas.) 


551 ± 48 


410 ± 42 


145 ± 25 


18 ±22 


528 ± 43 


MC(pred.) 


489 


362 


126 


33 


456 


N(S*-)/N(A) 


uncorrected 


DATA (%) 


2.6 ±0.8 


2.2 ±0.9 


3.7 ± 1.5 


2.9 ± 1.0 


2.4 ± 1.1 


MC ( fi ) 


1.5 ± 0.1 


1.6 ± 0.5 


l.G ± 0.8 


0.1 ±0.5 


6.5 ± 0.5 



A summary of the number of extracted TP events and the E°/A ratio is given 
in Table 16. 



5.3 Yields of strange resonances and heavy hyperons 

The integral yields of strange resonances and heavy hyperons produced in 

CC are computed multiplying the results of the fits by the ratio MC(pred.) /MC(meas.) . 

The results are presented in Table 17 as fractions of produced by heavy 
strange particles and resonances. 
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Table 15 



Att summary 



N(S-) 


full sample 


A fragmentation region 


type of target nucleon 






XF < 


> 


up 


vn 


DATA 


42 ±30 


21 ±24 


18 ±17 


54± 18 


-11 ±24 




43± 18 


33± 15 


13 ±9 


9±8 


36± 16 


MC(pred.) 


60 


47 


15 


14 


47 


N(H-)/N(A) 


uncorrected 


DATA (%) 


0.5 ±0.4 


0.4 ±0.4 


0.7 ±0.7 


1.6 ±0.5 


-0.2 ±0.5 


MC (%) 


0.4 ±0.2 


0.4 ±0.2 


0.4 ±0.3 


0.2 ±0.2 


0.4 ±0.2 



250 - 




A/fA^Y) (GeV/c^) M(A''y) (GeV/c^) 



Fig. 24. A7 invariant mass distributions for both MC (left) and data (right). The MC 

plot shows the expected signal peak and background contributions from — > Att^ 
and S*^ — > Att'' decays with only one reconstructed photon. 

5.4 Discussion 



The results of our study confirm discrepancies reported earlier [14] in the de- 
scription of the strange resonances and heavy hyperon production in neutrino 
interactions by the LUND model with default parameters [30,31]. These re- 
sults could be potentially used to tune the MC parameters responsible for 
the fragmentation into strange particles. Moreover, an additional analysis of 
events with multiple production of neutral strange particles could be very 
useful in this respect. Such an analysis is currently in progress. 
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M(n) (GeV/c^) M(Yi) (GeV/c^) 



Fig. 25. 77 invariant mass distributions for MC (left) and data (right). Both photons 

have been reconstructed as conversions in the DC fiducial volume and identified by 
our identification procedure. A clear peak corresponding to the vr" signal is visible 
in both distributions. The parameters Pq and P-j show the mass and the width of the 
Gaussian function after the fit. 

Table 16 



E A7 summary 



N(SO) 


full sample 


A fragmentation region 


type of target nucleon 








XF > 


vp 


un 


DATA 


29± 10 


17 ±9 


16 ±7 


16 ±7 


13 ±7 


MC(meas.) 


82± 12 


50 ±9 


37 ±8 


19 ±7 


61± 10 


MC(pred.) 


80 


57 


22 


22 


57 


N(SO)/N(A) 


uncorrected 


DATA (%) 


0.4 ±0.1 


0.3 ±0.2 


0.7 ±0.3 


0.5 ±0.2 


0.3 ±0.2 


MC (%) 


0.7 ±0.1 


0.6 ±0.1 


1.2 ±0.3 


0.5 ±0.2 


0.8 ±0.1 



Table 17 

Corrected fractions (in %) of observed Kg and A decays that originate from the 
decays of strange resonances and heavy hyperons in the NOMAD data compared to 

the dcfdiill MC prediciioiis. 





K*+ ^ K07r+ 




S*+ ^ A7r+ 


S*- ^ Att- 


^ A7 


S — Att 


DATA 

MC 


15.5 ±0.9 
31.4 


8.7 ±0.7 
13.1 


5.8 ±1.1 
16.6 


2.6 ±0.8 
3.9 


7.3 ± 2.4 
12.7 


1.9 ±1.7 
1.5 
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6 CONCLUSION 



We have reported the results of a study of strange particle production in i/^ 
CC interactions using the data from the NOMAD experiment. Our analysis 
is based on a sample of z/^ CC events containing 15074 identified K^, 8087 
identified A and 649 identified A decays. This sample represents at least 
a factor of 5 increase in statistics compared to previous (anti)neutrino exper- 
iments performed with bubble chambers. Yields of neutral strange particles 
(Kg, A, A) have been measured in this analysis as a function of kinematic 
variables. For A production such measurements are performed for the first 
time in a neutrino experiment. The decays of resonances and heavy hyperons 
with identified Kg and A in the final state have been analyzed. Clear signals 
corresponding to K*^, E*^, S~ and E° have been observed. This study is po- 
tentially interesting for the tuning of Monte Carlo simulation programs and is 
also of special importance for a quantitative theoretical interpretation of the 
A and A polarization measurements reported earlier [27,28]. 
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